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ABSTRACT 


Properties are developed for the failure rate function 
ang residual mixing probabilities resulting from a mixture 
of distributions. After a detailed investigation of these 
Dueperutes 10r purelywexponential mixtures, two more general 
mexEUres sare consaqered: Emphasis is placed on relating the 
current failure rate in the mixed population to the popula- 
maon's reetouedl composition. The results have potential 


application in determining an appropriate period of "burn in." 
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I. INTRODUCTION 


Sieeesce thab at time t=O there exists a population of 
Similar elements. This population could be made up of a type 
©f transistor, a type of hydraulic check valve, a certain type 
Dieplock, Cl cevera. But whatever the elements composing the 
population may be, they are all designed to perform the same 
task. And being a functional physical object, each element 
will have some random lifetime, depending on a known or unknown 
mrOototlty Gisturoupron, during which it. functions properly. 
When this random lifetime has elapsed, the element fails. 

The usual definition of the probabilistic failure rate 
PuMmecu lon Ol saiecuemeny With random lifetime will be used 
MieOuehOULetatowumesisS. Nemely, the failure rate function of 
an element is the density function for the lifetime of the 
element divided by the element's survival function. The 
scenario can now be further developed by requiring that the 
population consist of subgroups of elements, with each sub- 
group having its own unique failure rate function applicable 
wo cach element in tite subgroup. There is, then, a mixture 
of failure rates in the population. It will be assumed that 
at time t = 0 the exact proportion of the total population 
held by each subgroup is known either explicitly in the dis- 
Gmewe Case Ol amplicitly through a cumulative distribution 
mune rrOom in tMe continuous case. And so, the proportion of 
Uecwpeopulayion ay time t = 0 having a particular failure rate 


is known. 





Now that the Population has peeén defined, it is possible 
to ask the two questions whose answers this thesis is con- 
Camnea witha First, if an element is selected at random from 
the population at time t = 0, what is the nature of the 
mandem’y selected element's failure rate function? And given 
piesa particular element selected at random from the popula- 
Prem ab Cime t = O° 1s still functioning at time’ t > 0, what 
mS he prebabiligy thas that element was a member™ofr a 
Speecific subgroup? 

For someone wanting an explanation of the practical impor- 
lmace Of CThese questions, a plausible example might be the 
following. SsUuppese that the populaticon™consPsts “or transisvors 
which were not marked as to subgroup membership due to some 
overoucht Gurmasene production process. If the’ installation 
of one of the transistors into a piece of electronic equipment 
ae time t = 0 as required, What can bé™predicted about the 
Plime ravewor= ene installcdrtransistor? "And on the basis 
SL PCs eperiokmance, What can De™sald@about’thatewtransisver’s 
Subgroup source? 

Others have already developed some of the answer to the 
first question raised above about the nature of the failure 
rate of an element selected at random from the population at 
time t = 0. For instance, Barlow, Marshall, and Proschan 
(Ref. 1] have shown that a mixture of distributions each 
having a nonincreasing failure rate itself has a nonincreasing 
failure rate. Proschan [Ref. 2] applies this result to help 


explain an observed phenomenon in the aircraft industry and 





to conclude that a mixture of exponentials will have a 
eee a sine failure rate. Esary, Marshall, and Proschan 
[Ref. 3] show that a mixture of distributions with decreasing 
hazard rate average has a decreasing hazard rate average. 
And Gnedenko et al. [Ref. 4] indicate that in a mixture of 
exponentials the failure rate always approaches as a limit, 
wath increasing time, the least exponential parameter value 
of all elements composing the population. This thesis sheds 
further light on the nature of the failure rate of an element 
drawn at random from a population composed of elements with a> 
mixture of failure rates. No claim is made that the exposi- 
Poitier eubCriMegmeomeive a Complete answer, and areas needing 
further study are indicated. 

Phew secenemireoation, awour the conditional probability 
that a randomly selected element has come from a subgroup 
with a particular failure rate function given that the element 
Mecvlllmuneumening at a icestain time, does not seem to have 
been piven much attention in the literature. And once again, 
the results presented here are certainly not exhaustive. More 
research is required if the question is to be completely 
em@owered. 

Briefly, the procedure followed in this thesis is to show 
results for increasingly more complex mixtures of exponentials. 
Omee this Nas been done, some attention is given to more 


general mixtures. 





IT. THE CASE OF TWO EXPONENTIALS (DISCRETE PARAMETERS ) 


bemmene Deopulation at time t = 0 consist of a proportion 
Py of elements distributed exponentially with parameter hy 
end a2 proportion Do of elements distributed exponentially 


with parameter i The equality Py + P5 = 1 holds. Assume 


2° 


that >. is less than A~. It is a well known result that the 


il 2 
failure rate functions h, (t) and h,(t) for elements in the 
two respective population subgroups are given by h, (t) = 
Syauel h, (t ) = Xo (see Figure 1). 

At time t = 0 an element is selected at random from the 
population. The following definitions are made: 
A, is the event that the element selected comes from 
the subgroup with failure rate hi: Ay is the event 
that hc element selected comes from the subgroup 
mith failure rate do: T is the lifetime of the 
element selected at random from the population. And 
Bits fr(t), and h,(t) are the element's survival, 


demisdityseand faijluresrate functions, mespectively. 


Using these definitions, some equations can be derived: 


Pele b) = P(T>t{A,)P(A,) foe Cl tA) Pia 
bs —A t -rat. 
pie +, Pe : 
a ar = —Ajyt —roat 
fip(t) Bie (i) p,A4e + pyroe 





h(t) = d- 


peace) = 2X 


Two Exponential Failure Rate Functions 


Figure 1 





(Ay - Aa)t 





— fp Py Ay * PQADE 
hp (t) ae 7 GL Se DY 
Bp) Py * Po & 
Nae 
—p,p,(aA, - dr Me ec 
. 3 


[Pa opel 92) 


(0] -- ho )t 
Palen og 


be’ '(t) = p.p.(aA, —- Xx y3_01 A2)t SS 
uP | aa: 1 5 
(Az - Aa)t 
P, + Poe 


Using the definition of h(t), the following theorem 
Results. 
Mneorem 1. For 0 <p, <1, Eoew ice latwensoip A 


<hyp(O) <A hedles . 


u zg 


Proof. Note that 


Am (0) = PiAq + DoAz = (] = Po)Ay + Pyro = a + poQ, - 1): 


Mim ole  bheoremuic va Special Case of the general 
Pecult Ssmeown an Ret. 1. 
Theorem 2. The failure rate h,(t) is a decreasing function. 
Proof. From the expression for hn’ (t ) given above it is seen 


that hn’ (t) < 0 for every t. j: 


Reference 4 indicates the general result of which the 
following 1s a special case. 


Ggeerem 3. lim h-(t) = 


ee it 


186, 





Peaeet. Whe result follows dimectly from the expression for 


A(t ) and the assumption made above that Ay <Ad- f 


Theorem 3 is an interesting result since it means that, 
no matter how great the probability Po of choosing an element 
Seecimne + = 0 trom the subgroup with failure rate ho may have 
been, the failure rate An (t ) eventually becomes as close to 
hy as desired. An intuitive explanation of this phenomenon 
is that as the element chosen at random from the population 
is observed to have an ever increasing lifetime without failure 
it becomes ever more likely that it is an element with a low 
failure rate. But the lowest failure rate it can have is 5 
Since there were only two possible failure rates in the popu- 
rat LO ronmwinmen the element was taken, and aT was the least 
on ULnese two. 8 lauuatively, itmwould seem that the longer the 
element functions the more likely should its failure rate be 
ro: Mie wets precess OF Theorem 3 reflects this Mkelihood. 
The range of h(t ) is bounded. 


Theorem 4. <hp(t)<’, for every t. 


aT 2 
Poot. — liemprooL tollows from. combining the results of 
Theorems 1, 2, and 3. # 

Some general statements can be made about the shape of the 
faire rave curve. 


Theorem 5. For Py <Po> the function hp (t ) has an inimlec ion 


re ee rrr ccm mee eee 


jPa 





Proof. THe function Ap (t ) has an inflection point when its 


second derivative vanishes. 


r r p 
h,''(t) = 0@Sp, -p e | La = 20 = 0oSt = ———— ln = 
34 aa 2 Dp 
ho) 2 
For Py <p5, this is a feasible time. That is, it is not a 
meeavrvye dnc’. 

Since Py < Pos it is seen from the expression for Ap! ' (t) 
Meavevais tunction has a negative value at t = 0, and it is 
monotone increasing for t>0O. Therefore, A(t ) is concave 
Berore the time at-which its inflection point occurs, and it 
is convex afterwards. # 


The value of the failure rate function at the time the 


ef 
AT r 


2 


ime lection point occurs is given by Ay Ct » where 


ee 
p 

"2 Oy ; my '#) 

As a result of Theorem 5 the graph of hp (t ) for the general 
case of a mixture of two exponential failure rates has the 
form shown in Figure 2. An intuitive explanation for the 
emepe Ol Uris cumeve aan De civen.. Prior to the time at which 
the inflection point occurs h(t ) always has values closer to 
ho than AW: lt seemenunhav al Ghe element chosen at randem 
Miesea relavively brief lifetime it probably has a high 
Herve waee, “Ime Concavaicy Of the curve, its shape tendency 
towards or prior to the inflection point seems to be an 
indicator that the element is more likely to have come from 


wme Subgroup with failure rate A ney tf aides sear eek 


a 2 


ae 





hp (t) 


—_—_ena ae ew eel Cl 


oso ean ff SBP ist ee es we Se fF = 


Failure Rate Function Resulting From a 
Mixture of Two Exponential Types (p, < Po) 


Figure 2 


aR 





peems to have a stronger attraction for the curve than hy 
eee prvar tosthesanflection»s point. 

AG ene amr lection point, the failure rate function hyp(t) 
has been shown to have a value exactly midway between dy and 


r The intuitive explanation for the inflection point 


5° 
occurring when this happens is that the element is probably 
equally likely to have come from either subgroup. The 
failure rates 5 and ho seem to have an equal attraction for 
ae Curve. 

Rigger The timewune Inflection point occurs, the curve 
becomes convex. The intuitive explanation for this is 
mnmebar LO thavweeinven above for the concave portion of the 
curve. 


For the case when Po <P> ne!) (Ge) As Zemeoronly at 4 


a 
iepavive Viienmeae sO Loreal! feasible times, t > 0, the 
function hp (t ) is convex and closer to Ao in value thaweacio 
do: The intuitive explanation for this involves the same 
reasoning used above for the case P, < Po: The case Pj = 
ollie: is the same except that h(t ) has an inflection point 
av t = 0, and it then becomes convex. 

A precise understanding about the conditional probability 
that an element randomly selected has come from a particular 
one of the two subgroups given that the element is still 
fimecvioning at a certain time requires the definition of a 
Pesoluuat mr xine Drobabllity, called a "mixing probability" 


Hencemect.  tnough the definition in this section is for the 


Gugeweceewomecdvials case, it will in later sections have an 


14 





easy extension pomuere Conolex Waepures. By the mixing 
probability at time t is meant the function P(t), ie ores 


such that P,(t) = P(A,|T>t), t>0. This may be rewritten as 


ee ae) ) ECR ese AS EGA; ) p,e ie 
Py(t) = PCAIT>) © —sepsey— = Ot) ay 
F(t) 


mer Che particular case of i = 1, this means that 


—A7t 
P_(t) 7 ae 1 
t) = = 2 
il —Azt ene O 7 - ASHE 
Ie), + Poe Py + Poe 
familar ly , 
ae 
Po lenia= : 
o (r90 - apre 
Po +t pie 


Note that P, (0) eee Gl P,(0) =p 


dh ae 

Once the mixing probabilities P, (t) and P(t) have been 
Spuecceamidmune above forms, it is easy to see that the 
intuitive explanation given above for the shape of the hp (t ) 
cumye iS@e@errect. Since Ay < Ads P(t) is monotone increasing 
to the value 1, and P,(t) is monotone decreasing to a value 
of 0. And for the case where hp (t ) Meise an inf lectwon point 


(the case when P) <Po); it can be seen that the following 


DetarLeonecmanps hold: 


il 





Pc Po iyedt tb < ty 
P,(t) = P(t) if t = ty 
P, (t) > P(t ) Po ot. ty : 


the maxing probabilities can be used in another way to 


see how the curve h(t ) eho se 


—A 7G is Porne 2" 


hy (t ) a AP, (t) + AoP,(t) = ECA). 
-7t -rot 
e + Doe 


mere 1S 90gemeondom Variable taking values aT and do at 
time t with probabilities Pit) and P,(t), respectively. So 
emcee URese s@endatuonal probabilities are known, the function 
hy (t ) is determined as an expected value. This is another 
fact having extension to cases more complex than a mixture 
of two exponentials, as will be seen in later sections. 

The first derivative hn’ (t ) can be represented in terms 


of expected values at any time t too. 


2 } 2 
EDC US) Comoeeme se (g.e 1°) ( 5 A aomeme aD 
yd i z je. Jd 


dt 
( 
J 


“Ajtye 


Il ot Nor FR 


p.e 
iJ 


{E, (A) - A, TP, (t) Home. = a. 2. 
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And so 


hy’ (t) 


fas is 
The 


concise 


Dn’ ! CES) 


t = t 1 as 2 2 
E, (A) = iP Cu) AP, (t) = -{E, (A )- Ey Cry} 


= - Var, (A). 


another way of showing that hp’ (t) <0. 


second derivative Ay’ (t) also can be put in a 


expected value form. 


2 r(ac 
= -{E,'(A") - 2E, (AE, CA 


t 
2s 2 
= = J a d 
a Bet) te (t )} 
2 2 : 
5 * | 
_ P| ~ 2r,E, (A) TA, - E,(A)}P, (t) 


= 2 

= E, {LA = 2hE, (A) ILA = E, (A) J} 
_ 5 2 2 

= E, tA = Bh E, (A) 4 hE, (A)} 


- 3 
peer ECA) 


The expression of Ap''(t) as the third central moment of A 


Will also have extension to general mixtures of exponentials. 


Ie 





ii, ep eeor THREE XPONENDTALSIYPES (DESCRETE PARAMETERS) 


Powo@mecemotdcrmne in detail the general case of a mix- 
pare Of Nn eCxpemential types, some results which have been 
fommad Lor a2 population composed of three exponential types 
are presented. Function definitions for the n exponentials 
@msec are given first, and then in the remainder of this 
Section it is assumed that n= 3. The more general defini- 
tions are needed in the next section. The definition of Ass 
fae eine paraitems that given in the last section 
for A, and A,. Also it will be assumed that the failure 


a a 
mace TunctiLonsmer them=pepulationssubproups satisfy 


Ay Sk SS emeere ae 
n n ae 
Be@t) = P(T>t) = © PMT mt te) PO, ) = @ pee “1 
il a iL 3. . i 
{=1 yo 
n 
_ mt a , ~wmAGt 
fin(t) Bp (t) = 2 pads 
i=l 
n 
» p,r,e °° 
i=1 
Ap (t ) - = ek. 
YY pe 1 
i=1 


When the definition of the mixing probabilities P, (t) 


given in the last section is extended to populations where 


18 





Coe Member of subgroups is n>2, the failure rate function 


can once again be written as an expected value. 


n 

An (t ) —, A,P, (t) = nN) 
i=1 

Using this fact and changing the upper index on all summation 

ees in ene NasGseceton from 2 to n, it is apparent that in 


me general case 
' = 
ha'(t) = -Var, (A). 
ha''(t) = E,{A = E,(a)}? 
Jk t u . 


A program was written in FORTRAN and run on the IBM 360 
as a means of observing the behavior of hy (t ) during various 
eume antvervals. For n = 3, examples checked with the aid of 
the program led to a conjecture which has been proven to be 
Prue . 

Theorem 6. If the population is composed of a mixture of 


three exponentials, with A, = A, + a, d. = le: 2a (a>O0), 


and Py £P; then hyn (t ) Mesmen inilection point at 


al ee : 7 
ty oF 1n( 53) . And hn»(ty) = do 


Proof. Note that at time t which is feasible since Py aera 
asia 


ie? 
P,(t,) = P.(tr). 


Piiere Lore, 
h(t, ) = ee = (A. - a)P, (ty) + AoP, (ty) + (A, = a)P. (ty) 


Ay {P, (ty) + P,(ty) 1 P(t) a NG 


1,8 





And also, 


'? 5 


{(Ay - a) - Ag}PPy(t) + (Ay -Ag)°PQ(t) 


5 
+ {(A, +a)- Lot Pa(t) 


5 E 
a {P3(t) = P,(t)} = Q. 


Theorem 6 details the case for n = 3 which is analogous 
to the n = 2 case, where it was shown in the last section 
that the inflection point, if there is one, occurs at that 
time t when hp (t ) is exactly midway between A, and A,. 
Theorem 6 describes the case for n = 3 when the inflection 
Meant wall oceur exactly midway between hy and Ao: And at 
the time this inflection point oceurs, for both the n = 3 
and n = 2 cases, the least and greatest indexed mixing 
probabilities are equal. Even the forms for the two cases 
©6l the times®: which @his®occcurs are analogous. For n = 2 
it was shown that, for P, <Po> chic Thileet ron pore OCceuics 


when 


p 
t= bay (F). 
2 1 Py 
Theorem 6 says the inflection point occurs when 


P3 


t = ore a = z eh 
2a Py Or. - ds) Py 


20 








It would be nice to find that hp(t ) always has an inflec- 


poem peimt when it takes on a value midway between A. and An’ 


al! 
pepomuiomomresulGwdocs Not carry over to more complex mixtures, 
ed 4G8tails to hold even for n = 3 when not all of the 


conditions of Theorem 6 are satisfied. As an example, 


consider the case where 


aT = 1 Par = 2 
do = 2 Po = 2 
Xba a eels 
3 P3 


Here 2.5 is midway between A, and i but the inflection 


1 cig 
point was found to occur when .24<t< .25 and 2.5739 <h,(t) 
Se ee. 
Theorem 6 does not extend to the case n = 4, for instance, 


even when conditions analogous to those required in the 


Poceremeare = satvstied. As an example, 


= ee ol 
d+ = 2 Der a wwe 
os ae Pa = +3 
hy = y oe = vt 


The inflection point was found to occur when .4469<t < .4470 
and 2.50955 < h(t) 2 CLS 

A particular property applicable here of a general result 
proven in the next section is that Ps (t) approaches the value 


1 while P(t) and P(t) approach 0 as t becomes infinite. So 


Call 





iste acaet or Lhercase n = 2, the least indexed mixing 
preobapmerty hase! as a limit with increasing time, and all 


other mixing probabilities have 0 as a limit. 


ee 





i, “Gatienal MExtuns OF BXPONENTIAL TYPES (DISCRETE PARAMETERS) 


The general case in which the original population con- 
sists of n subgroups, with all elements in any subgroup 
having the same exponential failure rate differing from the 
failure rate associated with every other subgroup, has some 
interesting properties. Definitions of Pp and P,(t) Carry 
over directly to the general case Peom €be n= 2 and n= 3 
cases. In the last section expressions were given for eG), 


A(t ) Javed irs eee cl axel ht" (t) for this mixture of n exponential 


dh 
my Des. 
Theorem 3 is just a particular case of a property pos- 


sessed by a general mixture of exponentials. 


Maieorem 7. limh_(t) = A 


F->00 q ily 
EOoOr . 
n 
(AZ - Ai )t 
PyAq a a pase 
h(t) = # 
T es (Ay - 4) 
p, + = p,e 1 1 
1=2 


Reference 1 shows that the function hy (t ) is monotone 
decreasing. This of course means that An! (t) <0 Tonia. 
When n = 2, Ap (t ) has at most one inflection point. 

This result has not yet been proven for the general case, 
even though computer runs of examples with various values of 


n have all revealed at most one inflection point for h(t). 


2S 





Mmieotb could be proven, that hp''(t) is monotone increasing, 
then it would follow immediately that h(t ) has at most one 
iM—~toChmon DOlms. But this does not hold in general. 
Consider, for instance, an example where n = 5 and the param- 


eter values for the mixture are 


id 


aT = 2 Py = zak 
ho = 5 Py = seals 
A2 = y P32 = mil 
Ay = 5 Py = -3 
le = 6 Ps = cu 


Here An’ ' (t) is negative at t = 0, increases to a positive 
maximum at t = .824, and then begins to decrease. All that 
Gan be conciuded is that if h(t ) is concave at t = 0 it 
must have at least one inflection point because An (t ) 
approaches 4 as a limit as t becomes infinite. 

A study Of mixing probabilities in the general case 


reveals a number of interesting properties. Writing P,(t) 


in the form 


Py 

Ps (t) = n 
2 AO ae a 
jel 
it can be seen that 

iP siOr aie.) 
igen P.(t) = 
eo * Se for = see ia 
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This is an intuitively plausible result in light of the 
limiting behavior of hy (t ) discussed above. The longer the 
element selected at time t = 0 continues to function, the 


Cl@screaics fawiure rate An (t ) Ces) FO =X its Creasonable, 


1 
then, to expect that the conditional probability that the 
element selected came from the subgroup with failure rate 
Ay> given continuing survival, should become greater with 
increasing time. And since the hn(t ) curve, for large t 
ea least, continually moves farther from the failure rates 


o> ‘5 eee if is reasonable that the conditional proba- 


ae 
bilities that the element selected has come from a subgroup 
with a failure rate greater than 5 should become less and 
ess. 

Additional facts about the behavior of the mixing proba- 
Pei hes Cam oe det cri ned. 


Theorem 8. If h,(0)>A,, P,(t) reaches a maximum when 


hp (t ) =), . 


Proof. An immediate extension to the general case of a 


result shown for n = 2 is 


dP, (t) 


ian aa. ial 


This means that 


dP, (t) 
dP, (t) 
ae 5 ee es 


ee 





dP, (t) 


=n 0 > An (t ) < As if 


The last line of the above proof shows that if hp (0) < ds 
tiie [UMC ha On P(t) has its maximum at t = 0 and is monotone | 
eo eCreacine Waereatter since h(t ) is monotone decreasing. 

The intuitive explanation for Theorem 8 seems to follow 
from the monotone property of hp(t). Since hp (t ) is always 
decreasing, it must, if hp(0) >A, ,5 get closer and closer to 
assuming the value hs until the time when Ap (t ) = Na: iGe 
Comoe CxPeCctTednewumen, Senat the cemditional probability the 
element chosen came from the subgroup with failure rate hes 
given survival until time t, should be increasing to a maxi- 
mum at the point in time when Nn (t ) = as As time continues 
beyond this, h,(t) gets farther and farther from A. It 
would seem it must be getting less likely that the element 
Smescn came Irem the subgroup whose elements have failure 
oe hae 

A property expressed in the proof of Theorem 6 for the 
cease n = 3 does not carry over to an arbitrary exponential 
mix. In Theorem 6 it was shown that ee LS aecascu nein 
n = 3 such that when A(t ) takes on a value midway between 
hy and dn UHe Cova l mixineeperobabalitysmass for ds above 
iio veal we IS eedha FCO the total mixing probability “mass 
below the midvalue. That is to say, when A(t ) reaches the 
midvalue, the element is equally likely to have come from a 


subgroup with failure rate greater than the midvalue as it 
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Pomcow@ave Come Irom one with failure rate less than the 
midvalue. But the example for n = 5 used earlier in this 
Seceuwon shows this result does not carry over to the general 
case, even though the A'S form an arithmetic progression 
just as the i,'s in Theorem 6 did. In the example, hyp (t ) = 
3 at t = .38276. But mixing probabilities at this time for 


rs values below the midvalue i’. satisfy P,(t) + P,(t) = ,42433, 


3 


while mixing probabilities for hs VvesooveCmnac Mlcwmarl uc 


satisfy 2) 2) ie P(t) = ,45834., 
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V. CONTINUOUS PARAMETER EXPONENTIAL MIXTURES 


Sem, tie bly Populatvwons which have been considered 
Gremumeoemurtn 2a Gascreve number of subgroups and, thus, a 
discrete number of hs parameters. In this section it will 
again be assumed that each element in the population has 
exponential lifetime indexed by a failure rate parameter, 
but now the parameter will not be peseeietednes discrete 
Votes woe ee iiaread, 1b will be assumed that it 1S a con- 
tinuous, nonnegative random variable A with cumulative 
distribution function G OX) and density BOA). When GOA) 
14s a step function, with G,Q,) = 1 for finite A» the result 
is one of the cases already considered. 

Vo VikeeOnommin viiS COMVIMNUOUS Case are extensions of 
tie CMesmiememoerete Cases. Hor instance, 


F(t ) —@iee >) = f P(T>tlaA = iene eda) we 
a a 


—~ rt: 
: g,OA)dr, 


where a>O and G, (a) = 0. 

A number of the general properties of Ap (t ) are the same 
here as when the parameters were discrete. Reference 1 shows 
thas h(t ) is monotone decreasing. Reference 4 indicates 
that =n h(t) is the minimum value that the random variable 


A can assume. And at time t = 0, h(t ) is the expected 


value of the random variable A. 
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Theorem 9. hp (0) = E(A). 


Proot. Nete that 


Cc 


F(t) = f e Fag (a) eee 


a 
= m™ me —-AtT 
fn(t ) Set (t) = E(Ae ) 
meen) 
ha(t) = ———— 


And therefore 
hp(0) = E(A) . # 


The properties of ny (t ) are illustrated in the following 
example. 
Example. Suppose that A is uniformly distributed over the 


interwed Ta, bi]; ae 0. Them 





= Nf 
F(t) = - emt oi “ eat _ Q7bt 
Ak : b-a ebe —eoeler 


(1 + ab yer 2! =i + fae” 


ace) = 
(b Se ea 


(1 + aioeme ae Ct yee OF 


J Ge a 
Ab (eat _ hs 


ag 





Using elementary limiting arguments, it is found that 


hn(O) = lim hp,(t) = 3 ee 


£0 c 





And also 


dm nn (t ) = a 


C0 

The analog in the continuous case of mixing probabilities 
1s 42 comditiemal density fumetion, g,0A). AM oils) albbelke akvoyal 
Pivyes the conditional density of the random variable A, 
given that the element selected at time t = 0 is still 
DUMC honing sotatiaame. F. wor he and heey Such that 
AA = Xs ae Ay is small, {ge (A) FAX will, then, give a good 
approximation to Che probability that an element surviving 
ag Game vs one with failure rate A satisfying 


Re SN ieercOnditwomal density 2s cdetanedspy 


ale i+l1° 


=\t 
e g, 6) 


oo 


g, Cr) = 


fe ae (u)du 
a A 


Heeememat, aS 1n Giserete cases, the failure rate func— 


tion can be interpreted as an expected value at any time t 


Since 
f re aa (2) 

hy»tt) = Bo dg, CA)GA = E(A|T >t) = E, (A). 
f e Fag (a) a 


a 


2) 





Theorem 8, too, has its continuous analog. 
Micorenw). ht ACO) > 2, g, 0) reaches a maximum when 
h(t) = A. 


Proof. Note that 


co 


Pores Cs ee (u)dul[-re7 "ge (a) ] + fom Onis Wen eCaul 
i ase a 
* 2 
7 ty ep (u)aul 
a 
= fhp(t) - Ade, Qa) . # 


Angsot commese, as before, if hn(O) <r the function g, 0) 
has its maximum at t = 0 and is monotone decreasing there- 


after since nn (t ) is monotone decreasing. 
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VI. SOME NONEXPONENTIAL MIXTURES 


in this section the requirement that all elements in 
the population have exponential failure rate is relaxed. 
Although a completely general mixture is not considered, 
previous theory is extended to mixtures more complex than 
the exponential ones considered thus far. 

A useful function in studying more complex mixtures is 
the hazard function. This function has been used by others 
(see, for instance, Ref. 3). The hazard function of an 
element with random lifetime is the function H(t) defined 
by H(t) = -In F(t), where F(t) is, as usual, the element's 
Sumvivelammevton.  Andait f(t) and h(t) are the element's 
Gensivy and tTaiiure ware functions, respectively, the follow- 


ing relatuemsnips hold. 


n(t) = £4) 2 y(t) 
ECt) 
Map er (tL) = Haye? 
And since 
F 


Hee f hesods , 
0 


it £follvems that 


Ae = nteyer"*) 


Note that H(t) is an increasing function since F(t) is 
Geecmeocine, 
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Now suppose that the population is composed of n subgroups. 
All elements in subgroup i have density function f(t), 
survival function F,(t), failure rate function h, (t), and 
hazard aunc tion H, (t). None of these functions for a sub- 
Broupeus Hqual Co fts corresponding function in any other 
Subgroup. And again at time t = O each Subgroup i makes up 
e propore.on Ds of the total population. Without any other 
Demmi culenmc. cas Would describe the most general mixture 
of a discrete number of subgroups. 

Consider now a case where each failure rate h(t), sles Le 


is greater than h, (t) burawcomovanl avy alll tames t. “hat aise 


h, (t) = h, (t) ge 


AL 

h,(t) = h, (t) +. kK, 

h,, (t) = h, (t) + a 
where 0 = ky < Ky So gee Kae Then 

19 

H, (t) = ‘ h, (s)ds = H, (t) + k.t 
And since h, (t) < h,(t) ie ae h(t), H, (t) < H,(t)< ...<H(t). 
And also 


~[H, (t) + k:t] 


f(t) = ny(t)e Mt) = Ch (t) + k, Je 


ame) eee) = eeLHy (tb) + kyt) 
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An elemeny selected at random from the population at 
time t = 0 will, then, have random lifetime described by the 


failure rate function 


n n 
5 pm.) Peoria (t) + k. Je7lH1(t) + kyt] 
asl Se th 4 
hyp (t ) = i ee oes) ee ee 
n n 
Oe , (6) me LEAT D + k:t] 
i=1 eae 


And since F, (0) = 0 by assumption, it follows that 


An (0) = p,h, (0). 


ii 


Iu Ss 


al 


This corresponds to the result shown for exponential cases 
Ghat Np (0) is the expected value, weighted by the p,, of 
fairest vescmilme une population at time t = 0. In Pact, just 
as before, A(t ) can be interpreted at all times t as an 


expected value since mow the mixing probabilities are given 


by 
p.F. (t) Dp a walgeGte) Te kit] 
Baa) = — ie J. 
Sei ct ) z ee eure) Shae 
a J J spa J 
And so 
n 
An (t ) = 2 h, (t)P, (t) 


i=1 
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P,(t) Caonmackwualilvy= be put into the same form that the 


ferns probabality had in the exponential case since 


ees” 


fii s@@eens that the-limiting result shown for exponential 


CasSCaommewas here. “That 1s, 


I 
= 


itor 2. 
lim P(t) = 
1 nee = (0) oes al 


HI 
NM 
we 
we 
2) 


Another limiting result holds too. If hp (t ) Ls pueean 
the fezmn 
-k3t 


~K.e 
eae 


48, 


re Le 
A(t ) = h, (t) + aa 
ae 
i 


eS Sl SS 
NM 


HY 


p 


Jute 


it can be seen that = [h(t ) - h,(t)J = 0. Just as in the 
expcnential case, Am (t ) ultimately approaches as a limit the 
least tambure raze Lfunctien present 1n the original population. 

Theorem 8, too, has its counterpart. 


Theorem ll. If h,(0) > h,(0), P,(t) reaches a maximum when 


h(t) = h,(t). Otherwise, P,(t) is continually decreasing. 


3D 





Proof. The argument exactly parallels that used before, 


since 


dP, (t) 


—4— = {hg(t) - h, (t)3P, (t) 


And A(t ) SeGs continually closer to h, (t) since 


is a decreasing function. # 


Comsider now @ pOpulauaon Ii whaen ghe failure rates for 
Subgroups Satisfy 

h, (t) = k,h, (t) 

h(t) = koh, (t) 


where 1 = k_.mw« k. <... < KS Note that in this case 
Fi(t) = (F(e) 74. 

The result in the exponential case that hy (t ) is a 
decreacsms function cannot be expected to carry over at all 
times Here since h, (t) may at some times be increasing. 
However, hyp (t ) can be expressed as a time dependent multiple 
KCu stone h(t), and y(t) is a decreasing function. To see 


this it is helpful to consider the case where the multipliers 
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are not discrete k,'S but instead come from a continuous 
random variable K > 1 with cumulative distribution function 
G(k). When G,(k) is a step function, and Gy (k_) = 1 for 


some finite kK? the discrete case is described. 


Cc 


= 7 = k 
F(t ) = : [F, (t)] dG, (k) 

CT Goi laa = alludes (ee Pe dices) 

T : it al K 

hn (t ) = h, (t)y(t) » where the function y(t) is 


foe) _ in k 
lee Le Ge eC Uk = puowk 
1 1‘ ) 4 x ) ‘ ELKF, (t) ] 
W006 | se Sees 


= K 
s CF (t) Tac, (c) BER {tI 
1 


Note that 


coe) 


t h k e 
y(t) (CF) (t) Mae (ie) 3 


co 


= k Digs k-1 
- ; [F, (t¢) Jae, (k) ; KEP, (t)]> “[-f, (t) Jag, (k) 


Loe) oe) 


f KEK, (t) J“ae,,(k) f «fF, (6) FP E-e, (tb) JaG, (ie). 
1 


iiweemeams that y'(t) < 0 if and only if 


oo 


Cy KLFy(t) Vac, (k))© < ¢ (F(t) Jaa (ke) sk E F(t) )*ae,(k). 
a 1 


3H 





But 


eo] 


= 2 ~ = e 
U KUF, (t)Tac,(k)) = tr (WER (6) 7%) (xeVEF, (t) ac, (i) 
il: 


oo Lo @) 


= k ee k 
< / LF, (t)] dG, (k) : k CF, (t)] dG, (k) 


by the Cauchy-Schwarz inequality for Stieltjes integrals. 
This proves that y(t) is a decreasing function. And this 
means that once the hyp (t ) curve crosses below an h, (t) curve, 
it never intersects that h, (t) curve again. This result, 
too, was found in the exponentials mixture case. 

Returning to the discrete mixture, h(t ) again approaches 
the least failure rate function present in the population 
as time increases. 
Theorem le. 


hyp (t) 
lim h (ey) ~ 7 ° 


CF > 0 


Eeeeres Wore the definition of h(t ) it can be seen that 


: | 
(ky - ky)H z(t) 
h(t ) ese a a 
lim Sim —— ee ee 
im 2 es ¢ (ky - ky )HL(t) ss 
Py eal p5e 
i=2 F 


Noting that 


P(t) = 
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it can be seen that hp(t ) has the usual expected value 
tivwempretation at any timnest. That is, 


ha(t) = £ h,(t)P,(t) = E,[h,(t)] 


L 
L= 


Ih 
Theorem 11 holds here aiso since once again 


dP, (t) 
aan h, (t)}P, (t) 


NWeree Gemearch is needed Co determine the existence of 
possible necessary and sufficient conditions which failure 
rate functions in the original population need to satisfy 
in order that properties of h(t ) and P,(t) seen for mixtures 


cOns@@ered in this thesis hold. 


39 





BIBLIOGRAPHY 


Cogie pce ctomai. A. W., and Proschan, F., 
"Properties of Probability Distributions With Monotone 
Hazard Rate," The Annals of Mathematical Statistics, 
v. 34, p. 375-389, June 1963. 


Proschan, F., "Theoretical Explanation of Observed 
Decreasing Failure Rate," Technometrics, v. 5, p. 375- 
383, August 1963. 


Booty, JoePacse Marshall, #A. W., and Proschan, F., "Some 
homme miakiy Applications of the Hazard Transform," SIAM 
Journal of Applied Mathematics, v. 18, p. 849-860, 
June 1970. 





Gwet@em-One>. V., Belyayev, Y. K., and Solovyev, A. D., 


Meathemetical Methods of Reliability Theory, Academic 
Paess, 1909. . 


40 





ILO} 


eee DIESTRIBUTION LIST 
No. Copies 


Defense Documentation Center 2 
Cameron Station 
Alexandria, Virginia 22314 


babies Ode 0212 eC 
Naval Postgraduate School 
Monterey, California 93940 


NMepe@cmmbLiotressor J. D. BSary, Code 55 Ey i 

Department of Operations Research and 
Administrative Sciences 

Naval Postgraduate School 

Monterey, California 93940 


LT Michael Wayne O'Bar, USN al 
3220 Apache Avenue 
San Diego, California 92117 


Faoteomor Sam C. Saunders 1 
Department of Mathematics 

Washington State University 

Pullman, Washington 99163 


Professor A. W. Marshall i 
Depamtment of Statistics 

University of Rochester 

Rochester, New York 14627 


Professor Z. W. Birnbaum 1 
Department of Mathematics 

University of Washington 

Seattle, Washington 98105 


Peoressor R. FE. Barlow al 
Department of Industrial Engineering 
and Operations Research 
University of California 
Berkeley, California 94720 


Professor Ernest M. Scheuer al 
Management Science Department 

san Fernando State College 

Nevitnridge, California 91320 


Professor Frank Proschan al 
Department of Statistics 

Lie fiteerda State University 

Tallahassee, Florida 32306 


47 





ag! igs 


ae 


Chief of Naval Personnel 
Pers l1b 

Department of the Navy 
Yashemeuen. D.C. 20370 


Naval Postgraduate School 

Department of Operations Research 
avon iowrdbave oC Temces 

Monterey, California 93940 


42 





Secunty Classification 









DOCUMENT CONTROL DATA-R&D 


entered when the overall report Is classilied) 
2@.REPORT SECURITY CLASSIFICATION 


The Failure Rate and Residual Mixing Probabilities 
Resulting From a Mixture of Distributions 


(Security classification of title, body of abstrect and indexing annotation must be 






ORIGINATING ACTIVITY (Corporate author) 


Naval Postgraduate School 
Monterey, California 93940 









_- REPORT TITLE 











1. OESCRIPTIVE NOTES (Type of report and inclusive dates) 


Master's Thesis: September 


. AUTHOR(S) (First name, middle initial, laet name) 






Michael Wayne O'Bar 


. REPORT OATE 7a. ee OF PAGES 76. NO. ae REFS 
september 1972 


ja. CONTRACT OR GRANT NO. 98. ORIGINATOR'S REPORT NUMBER(S) 









6b. PROJECT NO. 


ec. 96. OTHER REPORT NOUS) (Any other numbere that may be aeaigned 
thie report) 


4. 
10. OISTRIGUTION STATEMENT 


Approved for public release; distribution unlimited. 


(1. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 


Monterey, California 93940 





13. ABSTRACT 


Properties are developed for the failure rate function and 
residual mixing probabilities resulting from a mixture of distri- 
butions. After a detailed investigation of these properties for 
purely exponential mixtures, two more general mixtures are 
considered. Emphasis is placed on relating the current failure rate 
in the mixed population to the population's residual COMP ost Gems 


phe results have potential application in determining an appropriate 
mr10d of “burn in." 





S/N 0101-807-6811 43 Security Classification... 


A-31408 








Security Classification 


Fr LINK A LINK £8 ; 
KEY WOROS 


Reliability 
Failure Rate 


Peaxture of Life Distributions 





D D HORM 1473 (Back) 


S/N 0101-807-68 4y 
a Security Classification A+31409 














The failure rate and residual mixing pro 


ina 





